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We use the Planck LFI 70GHz data to further probe point source detection technique 
in the sky maps of the cosmic microwave background (CMB) radiation. The method 
developed by Tegmark et al. for foreground reduced maps and the Kolmogorov parameter 
as the descriptor are adopted for the analysis of Planck satellite CMB temperature 
data. Most of the detected points coincide with point sources already revealed by other 
methods. However, we have also found 9 source candidates for which still no counterparts 
are known. 


1. Introduction 


New data by Planck experiment provide more accurate information on the cos¬ 
mic m icrowave _ backgro und 


(CMB) (Ade et al. l2013ac Bouchet & on behalf of the Planck collaboration! 20141) . 
The latter is known to carry crucial information on the properties in the early Uni¬ 
verse, while the temperature fluctuations in the CMB maps are basically Gaussian 
( Komatsu et al. 200.4 Fergusson et, al. 2012 : Ade et al. 2013c). Along with that, 
certain non-Gaussian properties as the Cold Spot, the low multipole anomaly, etc 
( Cruz et, al1l2004 VielvalEoiO ; de Oliveira-Costa et all 2004) have been als o noticed, 


including in the Planck data (iGnrzadvan et al.ll20l4TRassat et al.ll2014h . Among 
the important issues regarding CMB, is the separation of point sources (galaxies, 
quasars, blazars, GRBs, etc) in the data. Here we study this issue in the Planck 
data using the method of Kolmogorov stocliasticity parameter. 

The cosmic microwave background (CMB) maps have been examined for the de¬ 
tection of point sources as potential foregrounds, which coul d be either cosmol ogical 
or Milky Way objects emitting thermally or non-thermally ( Ade et al. 2013al b'l. A 
corresponding catalog is made available in NASA/IPAC Infrared Science Archive 
Various methods, including the wavelets and ne edlets, have been us e d for separation 
of point sources in pixelized sky maps (see, e.g., Curto et~ahl ( 2013 1: Scodeller et al 


1 


















































March 24, 2015 1:3 WSPC/INSTRUCTION FILE K-map39-mpla 


( 20121 ): Batista et alJ (2011); Lanz et ajJ (120131) ). Most of the detected sources in 
CMB maps coincide with know n radio sources, quasars, blazars, alth ough some 
sources remain still unidentified (jWright et al.ll2009t IJarosik et al.ll201ll) . 

The K olmogorov stochasticity parameter (KSP) has been alr eady involved for 
such aims ( Gurzadvan fe Kocharyan 20081 : Gurzadvan et al. 2010 ) using Wilkinson 
Microwave Anisotropy Probe (WMAP) data: some of the sources revealed by that 
method initially had no counterparts, however their counterparts have been detected 
later by the Fermi satellite as g a mm a so urces and were included in Fermi LAT 1- 
year Point Source Catalog (jAbdo et ah 12010 ). In the present study we continue 
the application of the KSP t echnique aiming to detect point sources in Planck LFI 
70GHz data set (jZacchei et al.ll201ll). and fo r that goal we borrow the maps cleaned 
by the method developed in Tegmark et_al. (|2003|). 

This paper is organized as follows. First we introduce briefly the Kolmogorov 
method, then we construct cleaned K-maps using the Tegmark et al. method, mod¬ 
ified for such maps. Then we use the resulting maps to reveal the point sources and 
other non-Gaussian structures on the CMB map. 


2. Kolmogorov method and CMB maps 


In his original work Kolmogorov Kolmogorov (119331 ) has introduced a method for 
determining whether a given random number seq uence Xj , i = L ..,7V, ordered in an 


increasing way, obeys a given statistic or not 


OOCj UOiiOO V - - X- ^ ^ J. T ^ 

jArnoldll2008allbU(]. l2009a. 


m . For such a 


purpose a so-called empirical cumulative distribution function Fn(x) is calculated, 


F n (x ) 


0, 

x < X\ 


k/N , 

Xi < x, k = 1,2,..., N — 1 

(1) 

1 , 

X N <x 



where k is the number of elements which obey to the relation Xi < x. Having 
assumed a particular theoretical cumulative distribution function (CDF) F(x), the 
parameter Aat is easily calculated, 

\n = VN sup \Fn(x) - F(x)\. (2) 


Kolmogorov proved that in the limit N —> oo, Ayr, which is a random variable, has 
a cumulative distribution function c f > (A) reading as 


+oo 

*w = E (-i) ie_2/A2 - ( 3 ) 

i=-oo 


The function 4>(A) can be expressed as a particular v alue of theta functions, since 


$(A) =i? 4 (0,e- 2A2 ) =t? 3 (0, 

— e 2x2 ) (see Abramowitz & Stegun 

1970. 16.27). 

In the series of papers 

Arnold 

2008a 1 

u 

2009a 

b|) Arnold 

examines the ap 


plication of this statistic to determine degree of randomness of a given sequence. 
Although this statistic was intended for i.i.d sequences, he applied this statistic 
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Fig. 1. Planck LFI 70GHz Cosmic microwave background temperature map with resolution pa¬ 
rameter n s ;d e = 2048 (top) and K-maps with resolution parameters n si< j e = 128 (bottom). It can 
be seen a round spot on the CMB temperature map affected by inhomogeneous coverage of the 
sky. This spot a ppears on K-map as a region with unusual high value of KSP (see Fig. [2 and 
Figs. 10,15,18 of lZacchei et all j201 il l 1. 

to arithmetical and geometrical progressions and revealed the relative degree of 
randomness e.g. of arithmetical progression with respect to geometrical one. 


3. Cleaned CMB K-map 

3.1. Planck LFI 70GHz CMB and Kolmogorov maps 


It is well known that the CMB sky has basically Ga ussian distrib ution, al thoug h 
along with certain non-Gaussian features (see e.g. Bernui fc Reboucail ( 20 id ): 
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Ade et al. (I2013dj) : Komatsu & SpergeJL (120010 1. Based on this empirical fact, the 


Gaussian distribution is used in Kolmo gorov method to co nstruct a K-map. In the 
case of Planck LFI 70 GHz CMB map (IZacchei et ah 2011) with resolution param¬ 
eter n s ide = 2048, for every compact region of the CMB map containing 256 neigh¬ 
boring temperature pixels, one corresponding pixel value of the K-map is obtained. 
But for Planck LFI 70 GHz CMB map with resolution parameter n s id e = 1024, 
the same procedure is done for 64 neighboring temperature pixels. Therefore for 
both CMB temperature maps with resolution parameters n s id e = 2048,1024 we get 
K-maps with resolution parameter ?r s ide = 128. This adds some statistical inhomo¬ 
geneity in our calculations but provides more maps for Tegmark et al. foreground 
reducing method. For both Planck LFI 70 GHz CMB map s with n s id R = 2048 
and n s ide = 1024 resolution parameters (details in G orsk i et al. 02005 )). we ob¬ 


tain a K-map with n s id e = 128 parameter. This means that if the CMB map has 
»*ma P = 12ng ide = 50331648 pixels, then K-map has n map = 196608 pixels. This 
is due to the fact that KSP is a statistical parameter. Then the KSP distribution 
maps over the whole sky can be obtained for Planck LFI 70 GHz data. It is seen 
in Fig. |T| that the Galactic disk region has higher and saturated KSP values, which 
indicate that it has a non-Gaussian distribution, as distinct from the CMB. Also a 
lot of pixels have a high value of KSP, but most of them, as it will be shown below, 
are due to instrumental and other types of noise, also of non-Gaussian nature. 

Throughout the analysis we mainly use the HEALPIX program package for 
calculation of the spherical harmonic coefficients and to reconstruct the maps 


Gorski et al, ( 20051 ). For some tedious manipulating procedures with the spherical 


h arm oni c coefficients aim, (see Eq. JI]) below) we use the GLESP program package 
Doroshkevich et al. ( 20051 ). For maps with low resolution (in our case for K-map, 


n-side = 128,n map = 196608) the ai m calculation and map reconstruction the use of 
HEALPIX program package led to accuracy problems. To test whether the calcula¬ 
tion error is small or not for a low resolution map, we construct a unity HEALPIX 
map with n s ide = 128, and another one with n s ;d e = 2048, and run them through the 
aim calculation and the map reconstruction procedures. For both cases this proce¬ 
dure adds some non-isotropic noise which has almost zero mean (e) and a very small 
standard deviation (a). We obtain (e 204 s) = 10 -7 , (ei 2 s) = 10 -4 , (72048 = 5.6-10 - 5 , 
(7i28 = 9-10 -4 . Although some pixels around the poles have bigger error, the fraction 
of these pixels on the sky for n s id e = 128 is less than 0.5%. So for n s id e = 128, the 
error and standard deviation are sufficiently small, which enables one to calculate 
the aim and to construct the cleaned map for KSP using Tegmark et al. method. 
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Fig. 2. Planck LFI 70GHz foreground reduced K-map iTegmark et al . 1 <2003!) with resolution 
parameters n s ^ e = 128 (top) and regions with higher value of KSP (<f> > 0.5). A subtle star-like 
shape of inhomogeneous sky coverage pixels can be seen in both foreground reduced K-m ap and 
among higher KSP values of KSP indicated as yellow dots (see Fig. 5 of lAde et alj (l2013alP . Blue 
dots of high KSP values indicate those non-Gaussian regions that have counterpart in Planck 
Compact Point Source Catalog. Red (pointed by arrows) and yellow dots are those having no 
counterpart and regions affected by inhomogeneous sky coverage, respectively. 
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3.2. Modified Tegmark et al. method 

Any full skjy ma p ca n be_ represented via a series of Legendre spherical functions 
( Hivon et al.l 2002 ). where 


T{6, ^) = a lm Y lm{0 , ip), 
l,m 

aim = j T(9, (p)Y* rn (9,ip)sin0d9dip. 


( 4 ) 


As easily shown, the coefficients Ci of the Legendre polynomials Pfi cos 9) in the two 
point correlation function C ( 9 ) of the power spectrum are related to the a; m by 


C* — ( a lm a lm)i 

C(9) = -^^(2Z + l)C i P i (cos0). 


( 5 ) 


l,m 


So called cross power spectra is a common technique described in IHinshaw et al 
(12003 1. It can calculated via taking the cross-correlation power spectrum coefficients 
for different type of cq m , as 


C? = {a&a’J- 


( 6 ) 


Tegmark et al. ([Te g mark et al. 120031) CMB foreground reducing method uses cross¬ 
correlation power spectrum matrix to obtain weights for every band. This method 
is commonly used for obtaining a CMB foreground-reduced map from the original 
WMAP and Planck (SMICA composite map) CMB maps. In this method every map 
from different bands is weighted with weights w\. Here l is the multipole number 
and i the band index. It differs from the interlinear combination (ILC) method of 
weighting different bands suggested by the WMAP team (Jarosik et ah [20111 1 by 
dependence of weights on the multipole num bers. 


We use Tegmark e t al. method (ITegmark et al.1120031 : I Saha et al.1 l200fil 12008 ; 


Ourzadvan et al. 20091) to develop a cleaned CMB K-map from three couples of 
Planck LFI 70 GHz maps: nominal, nominal_ringhalfl, nominal_ringhalf2. Each el¬ 
ement in couple of maps contains both maps with resolution parameter n s ;de = 
1024, 2048. This method, based on the power spectrum comparison, assigns a weight 
w\ for each map i and multipole l. In order not to distort the original map, w\ should 
obey the relation 


2>( = L 


(7) 


But in our case it is mandatory to preserve original K-maps mean value because 
Kolmogorov parameter varies from zero to unity. Because first harmonic coefficient 
aoo depends only on mean value of map, we take aoo of cleaned map as a mean 
value over original maps aoo values 


aoo — (aoo) ■ 


(8) 
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Fig. 3. Angular distance between the nearest point source of Planck LFI 70 GHz compact point 
source catalog and KSP high value pixel ( d6 ) vs. Kolmogorov parameter <f> outside the Galactic 
Region |6| < 20. Red and blue dots are KSP high value pixels whose distance smaller and larger 
than 0.45, one pixel angular diameter of HEALPIX map with resolution parameter n s id e = 128, 
respectively. Black cross (+) over a red dot indicate K-map pixel around which too few observations 
are available. This pixel located over the spot is mentioned in Fig. [I] 

In the Kolmogorov method, if a sequence of random numbers T n obeys a theoret- 
ical distribution function F(x) then N realizations of this sequence gives <J>at, Aat. 
The remarkable point of the method is that, <f>jv has a uniform distribution and 
the mean value (<f>) = 0.5. Therefore for a CMB K-map the value 0.5 is a natural 
threshold for distinguishing non-Gaussian areas from others, since Gaussian tem¬ 
perature pixels cannot have $ > 0.5. Further, the cleaned K-map mean value and 
sigma are respectively ($) = 0.17, cr$ = 0.09, so pixels with $ > 0.5 mostly exceed 
the 3-cr region. 


4. Non-Gaussian areas and point sources 

4.1. Inhomogeneous coverage of the sky and cleaned K-map 


It is well-known tha t Planc k satel lite covered the sky inhomogeneously. It can be 
seen in Fig. 5 oflAde et_ah (|2013aTI that certain regions on the sky observed about 
thousand times but there are also sky regions that were not covered (observed) at 
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all. Kolmogorov parameter is a statistical parameter therefore it is sensible to such 
inhomogeneous coverage, too. So if a compact region of 256 temperature neighboring 
pixels is located on the border where number of observations (NOBS) decrease or 
increase sharply, then we have high value of KSP. These NOBS affected regions 
are colored yellow to indicate them among other regions with high value of KSP 
(<I> > 0.5) in Fig. [2j In contrary with these yellow colored pixels where number 
of observations decreases form very high values to mean value, there is another 
pixel of K-map around which hardly any observation was done (coordinates on the 
K-rnap l = 173.632, b = —50.48, colored red). So Kolmogorov parameter detects 
inhomogeneous coverage of the sky,too. 


4.2. Galactic region and point sources 

In Fig. [2] one can see the Galactic disk, the Large Magellanic Cloud (LMC) and 
other possible point sources. There are in total 1321 non-Gaussian pixels among 
196608 (about 0.7%). Only 153 among them are outside the Galactic disk region 
(|6| < 20), the Large Magellanic Cloud (l = 280.4136, b = 32.9310) and regions 
of sky affected by inhomogeneous sky coverage (so called ’’NOBS”). So, these 153 
regions are point source candidates. Indeed, most of them are found in the catalog of 


Ade et ah ([2011! ) (see Fig. [3] for more details). Note that even if one takes Galactic 


region as |6| < 18, then it would just add more 3 point sources to these 153. Hence 
these 153 point source candidates depend very slightly on Galactic disk mask choice. 
Among that 153 point source candidates only nine possible point sources remain 
unidentified in any of the known catalogs. 

We are comparing and confirming that our possible point source candidates 
have counterparts in Planck Compact Source Catalog in case when angular distance 
between a point source is equal or less than the average angular size of HEALPIX 
pixel. For a map of resolution parameter n^e = 128 the HEALPIX pixels angular 
size is d P i X = 27.5. But this restriction is very strict and some cases can be when 
point sources have large angular size. Therefore five of these unidentified point 
sources being close neighbors (69 < 1) might be also affected by Planck point 
sources. Remaining four unidentified point source candidates have counterparts in 
Two Micron All Sky Survey (2MASS) galaxy catalog. 
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_b 
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1.76 
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184.195 

45.784 
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0.49 

_b 

50.0562 

-57.0185 

0.504 

1.47 

14.676“ 

86.1328 

-38.6822 

0.95 

0.51 

_b 

59.0367 

-49.3128 

0.57 

0.47 

_b 

217 

-73.4963 

0.55 

4.66 

12.207“ 

237.682 

-48.9228 

0.54 

0.47 

_b 

a nearby galaxy exists 

in 2MASS 

catalog within radius of one 


arc minute 

b unidentified point source but very near to some object in 
Planck Point Sources catalog I Ade et al.l d2011f) 


5. Conclusions 

We use the Kolmogorov stochasticity parameter to probe non-Gaussianity of Planck 
LFI 70 GHz CMB maps. It is shown that for about 72% of the cleaned map, the 
Kolmogorov function has values within the interval 0.0 < $ < 0.2, which implies 
that the CMB maps are indeed Gaussian with high precision. 

Non-Gaussian pixels (with $ > 0.5) are rather rare (1321), i.e., less than 1%. 
Although this result is true for cleaned K-map, it shows that different sources add 
certain non-Gaussian noise into CMB map and they can be removed by foreground 
reducing methods. 

Another remarkable feature of both cleaned and original K-maps is that it shows 
some areas of inhomogeneous coverage of the sky by Planck satellite. It is obvious 
that some pixels in K-maps have high value of KSP due to this effect. 

For searching counterparts between these non-Gaussian regions and any avail¬ 
able catalogs, we utilize the main feature of HEALPIX map- equal area pixels. 
Therefore 1-pixel angular size is about 0.45 degrees in case of resolution parameter 
fiside = 128. It can be seen in Fig. [3] that this assumption works very well. K-map 
pixels outside Galactic disk region (\b\ > 20) with KSP value $ > 0.5 indicate 
possible point sources and other non-Gaussian regions. Indeed 143 of them have 
counterparts in Planck Compact Point Sources Catalog of LFI 70 GHz. But nine 
non-Gaussian regions still have no known counterpart. Note that, a non-Gaussianity 
appears around the region of CMB sky with extremly low number of observations. 
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